Quantization of the Taub cosmological model with extrinsic time 
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; The paper addresses the quantization of minisuperspace cosmolog- 

ical models, with application to the Taub Model. By desparametriz- 

^jy ing the model with an extrinsic time, a formalism is developed in or- 

der to define a conserved Schrodinger inner product in the space of 
.£h ! solutions of the Wheeler-De Witt equation. A quantum version of 

classical canonical transformations is introduced for connecting the so- 

ctf ■ lutions of the Wheeler-De Witt equation with the wave functions of 

the desparametrized system. Once this correspondence is established, 
boundary conditions on the space of solutions of the Wheeler-De Witt 
equation are found to select the physical subspace. The question of 
defining boundary conditions on the space of solutions of the Wheeler- 
De Witt equation without having reduced the system is examined. 
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I. INTRODUCTION 



General Relativity is an example of parametrized system, i.e. a system whose action is 
invariant under changes of the integrating parameter r ("reparametrization"), this invariance 
being a consequence of the covariance of the theory. This means that in General Relativity 
there is no privileged time variable. On the contrary, in the ordinary formulation of quantum 
mechanics there is a time parameter besides the true degrees of freedom, and the inner 
product remains conserved in the time evolution of the system. This difference between 
General Relativity and quantum mechanics, known as the problem of time is one of 

the main obstacles for finding a quantum theory of gravity. 

The evolution of a dynamical system is characterized by the way in which its dynamical 
variables evolve as a function of time. In this formulation time is a relevant physical pa- 
rameter clearly distinct from the dynamical variables. There is nevertheless an alternative 
formulation of dynamics (parametrized systems) in which time is mixed with the dynamical 
variables |4]]. A parametrized system can be obtained from an action 5* (q^,^) which is not 
invariant under reparametrizations by raising the time to the rank of a dynamical variable. 
Let us start with an action of the form 



By identifying q° = t, p Q = —h one can rewrite the integrand as Pidq % = Pi(dq l j 'dr)dr , 
i — 0, ,n. In this way the extended set of variables are left as functions of some physi- 
cally irrelevant parameter r. The set {q l ,Pi} can be independently varied provided that the 
definition of p t is incorporated to the action as a constraint H = p + = 0, so 

yielding the following action 



where N is the Lagrange multiplier whose variation assures that the constraint does hold. 
This action is invariant under reparametrizations r — > t + e(t). The time variable t satisfies 
the Poisson bracket 




n 
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{t,H} = l (3) 

Once the system has been parametrized, it can be reduced using any time variable 
provided that it satisfies ([$]). This kind of time variables are called global times ||. In 
order to generalize this restriction let us suppose that we know a globally well defined time 
variable t = i (g^pj)which satisfies 

{t,H}\ H=0 =f(q,p)>0 (4) 

The important fact is that / has a definite sign on the constraint surface (it could also 
be negative). In this case the variable t is a global time associated with the Hamiltonian 
H=f- 1 (q,p)H. 

The constraint H = could also be expressed in a set of variables in which the Hamil- 
tonian H has not the form H = po + h. In fact we can perform a canonical transformation 

{q\Pi} = {q° = t, Po = -h,q»,p,} - {Q\P,} 

where now the time is hidden among the rest of the variables. In other words, a constraint 
of the form H = po + h can be disguised by scaling it or by performing canonical transfor- 
mations. 

One of the main properties of the Hamiltonian formulation of General Relativity is that 
the Hamiltonian is constrained to be zero, making manifest that General Relativity is a 
parametrized system. In cases like this one, in which the theory is an already parametrized 
system, the invariance under reparametrizations means that there is no privileged time 
variable. To reduce the system means to select among the dynamical variables a proper 
global time, i.e., a variable which monotonically increases along any dynamical trajectory, 
to work as a physical clock. In this way we can express the evolution of the canonical 
variables as a function of this physical clock. The first step to reduce the system is thus 
to perform a canonical transformation in order to find a set of variables {q l ,Pi} where the 
variables go = t is a global time. The Hamilton equations are 
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dt 

d% 
dr 
dp^ 



Nf 



Nf 



-Nf 



dh 



dh 



The dynamics of the system is thus undetermined unless one fixes a gauge, i.e., unless 



One of the mains approximations for quantizing General Relativity begins by reformu- 
lating it under a Hamiltonian formulation (ADM formalism ||) . Within the framework of 
this canonical formalism or geometrodynamics it is supposed that the Lorentzian space-time 
manifolds M are diffeomorphic to R x S where S represents a collection of spacelike hyper- 
surfaces S parametrized by a real time parameter t (foliation). The Riemannian metric gij 
of one of these hypersurfaces £ play the role of the configuration variable. The analogous of 
the configuration space R n is the space of all the Riemannian metrics called superspace. 
The conjugate momentum 7T 4 - 7 is directly related with the way in which the hypersurface S 
is embedded in the manifold M, i.e., with the extrinsic curvature of the hypersurface S. The 
covariance of the theory under general coordinate transformations is reflected within this 
formalism in the presence of four constraints per each point of space-time. The so called 
Hamiltonian constraint assures the invariance of the theory under a changing of the folia- 
tion, while the momentum constraints assure the invariance under a change of the spatial 
coordinates used to represent the spatial geometry of each hypersurface. The states of the 
corresponding quantum theory \l/ [g^] are functionals of the spatial metric g^ which satisfies 
the quantum version of the classical constraints in accordance with the Dirac method. The 
quantum version of the momentum constraints implies that the wave function depends on 
the geometry 3 g of the hypersurface but not on the particular metric tensor g^ used to rep- 
resent it. The quantum version of the Hamiltonian constraint is the so called Wheeler-De 
Witt equation. 

Many of the tentatives for quantizing General Relativity began addressing the analogy 
between the Wheeler-De Witt equation and the Klein-Gordon equation. In fact both sys- 



one chooses a physical clock. Choosing the gauge r — t means choosing N (r) 
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terns have Hamiltonian constraints which are hyperbolic in the momenta. The constraint 
associated with the motion of a particle in a pseudo-Riemannian geometry has the form 

H particle = 9^ (/) Pi Pj ~ VT? = (5) 

The space of solutions of the Klein-Gordon equation can be turned into a Hilbert space with 
a positive definite inner product only if the background is stationary. In this case the Hilbert 
space of the physical states will be the subspace of positive norm, this being equivalent to 
consider just one of the sheets of the hyperbolic constraint surface. Choosing the coordinates 
in a way that = ( =>• g 00 = g^) and calling = —g^g^ we can write (|5]) in the form 



Hpartide = 9°° (poPO ~ l^PpVv ~ 900^ (6) 

= 9°° (po ~ \jY v PpPu + 9oom 2 ^ (p + \jY v PpPu + fl-oo^ 2 ) 

In addition it is necessary to find a temporal Killing vector of the supermetric which also 
should be a symmetry of the potential term (this property could be relaxed to a conformal 
Killing vector). In this case the proper time variable to reduce the system is the parameter of 
the Killing vector. Otherwise there would be pair creation. In order to build a good analogy 
with the relativistic particle is also necessary to have a positive definite potential term for 
playing the role of the mass term. If the potential is positive definite, the momentum po 
does not go to zero on the constraint surface H = 0. This means that the Poisson bracket 
{qo,H} = 2g 00 po has a definite sign on each sheet of the constraint surface. If we choose 



Po + ^^PyPv + goofn 2 = 0, the momentum p Q} and so {qo,!!}, will be negative on this 
sheet (provided that g 00 > 0). The other factor has then a definite sign on this sheet playing 
the role of the function / defined in (|j) . In this case / will be negative, being this the 
reason why {q Q = t, H} < 0. But t is still the variable which monotonically increases on 

any dynamical trajectory because \t,H\ = 1 where H = -, H The quantum 

physical states can be obtained by solving a Schrodinger equation with a positive definite 



operator h associated with the Hamiltonian of the reduced system yj^p^Pv + goo 171 " 2 - As 
it was said before, to fix the gauge t = r implies to choose iV = ~. The relation between 
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proper time and the time variable chosen to represent the hypersurfaces of simultaneity is 
dT = f~ l dt where T is the proper time, this being a consequence of the way in which the 
space-time interval is expressed in the ADM formalism. In geometrodynamics there is a 
conformal Killing vector of the supermetric J7| but this vector is not as well a symmetry of 
the potential term. Besides this potential term is the spatial curvature |J, which can be 
negative in some regions of the configuration space. Thus it is not possible to associate an 
operator h with the square root, as one effectively does for the relativist particle. 

There are thus two main approaches for achieving this quantization program. One possi- 
bility is to quantize the system without reducing it. The resulting Wheeler-De Witt equation 
is an hyperbolic equation while the Schrodinger equation associated with the reduced system 
is a parabolic one. The former has then more solutions than the latter, so being necessary to 
define boundary conditions in order to select the physical solutions. Besides it is not clear in 
this approach how to define a conserved inner product without having reduced the system, 
i.e., without knowing which variable plays the role of time. Another proposal is to perform 
a canonical transformation in order to find a Hamiltonian of the form (|5p with a positive 
definite potential term independent of the variable go- The formalism of the relativistic 
particle can then be applied using q as a proper time variable. The system can thus be 
quantized by means of the corresponding Schrodinger equation associated with one of the 
sheets of the constraint surface. The Hilbert space of the quantum states can be endowed 
with the natural inner product associated with the Schrodinger equation. This approach 
has the problem that different choices of global time variables can lead to different quantum 
theories. 

II. PROPOSED FORMALISM 

In this work we will address the quantization of minisupersapace cosmological models. 
The quantization program will be as follows. We will start with a Hamiltonian constraint 
such that none of the variables is a global time. We will suppose that it is possible to perform 
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a coordinate transformation so that a subsystem depending on just one pair of canonical 
variables {qi,Pi} is separated in the Hamiltonian constraint 



H q = h qi (qi,pi) - h q ^ (q li ,p li )n = 2,..,N 



(7) 



where the Hamiltonian h qi (qi,Pi) has the form 
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din V 



dqi 



Pl + V( qi ) 



(8) 



with V (qi) > and h q > 0. In order to find a global time, one should look for another 
canonical transformation {qi,Pi} — > {Qi,Pi} i — 1, ...,N such that 

Q 1 = t = t(q 1 ,p 1 ) (9) 
P\=Pt = [h qi (qi,Pi)}^ 

Qp, = q» 

(i = 2,...,N 

In this way we could separate an extrinsic time t, i.e., a global time variable which is a 
function of both original canonical coordinates and momentum. The generator function for 
the canonical transformation and the corresponding momenta are 



Pt 

Pqi 



-sinh (t) [V 
dF l 



(10) 



dt 
dF l 



-cosht \V{qi)Y 



1 • i r T w M-i dV 

--su>u[v( qi )] ._ 



-- S inht[l/( 9l )] 2 ^-^ 
2 agi 



so that 



dqi 



2 Pl + ^ (?i) = V (gi) sinh 2 t + V = V cosh 2 i = p 2 



The Hamiltonian in the new set of variables has the form 

H Q=P 2 t - K (fy>2V) 
7 



(11) 



Thus, when the Hamiltonian h qfl is positive definite and independent of time t one 

gets a constraint such that the analogy with the relativistic particle does hold, and one can 
quantize the reduced model by means of the parabolic Schrodinger equation associated with 
one of the sheets of the constraint surface. 

Once the system was reduced and quantized we want to find out what kind of boundary 
conditions should be imposed on the solutions of the Wheeler-De Witt equation expressed 
in the original set of variables. As was pointed out, the hyperbolic Wheeler-De Witt has 
twice the number of independent solutions than the parabolic Schrodinger equation. It is 
then necessary to impose proper boundary conditions for selecting the physical solutions. 
In order to do that we will follow the lines of work used in || . Knowing the classical 
canonical transformation for reducing the model, its analogue in the quantum level can be 
defined. In Ref. || the conditions for relating the wave functions corresponding to a pair 
of quantum-mechanical systems whose classical Hamiltonians are canonically equivalent are 
studied. If one has two arbitrary Hamiltonians related at the classical level by a canonical 
transformation corresponding to the generating function F\ (q, Q) , the main issue is to find 
out what kind of integral transforms can be defined in order to relate the wave functions 
corresponding to each quantum-mechanical system. Generalizing the Fourier transform, a 
relationship of the following kind is proposed 



where F (q, Q) is not in general the generating function F\ (q, Q) for the classical canonical 
transformation. In Ref. || it is shown, however, that this function coincides in fact with 
the generating function for the classical canonical transformation when the Hamiltonians 
operators satisfy the condition 



where some proper boundary conditions in the integration limits are also assumed. If the 
canonical transformation cannot be represented by means of a generating function of the first 




(12) 




(13) 
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kind, analogous integral transforms and conditions can be denned using the corresponding 
generating function. The inverse of the integral transform fll2|) is 

d 2 F(q,Q) 



/+oo 
dq 
-co 



dqdQ 



e -iF(q,Q)Q E ( ? ) 



(14) 



The canonical transformation defined by fllCf ) does satisfy the condition ( |i~3|) . The function 
F (q, Q) coincides then with the generating function F\ (q, Q) . Once defined this "canonical 
quantum transformation" the physical solutions of the Wheeler-De Witt equation can be 
found by transforming the solutions of the Schrodinger equation. Finally the question of 
defining proper boundary conditions for the solutions of the Wheeler-De Witt equation 
without knowing how to reduce the system will be addressed. 



III. APPLICATION TO THE TAUB MODEL 



A. Desparametrization 

We will study the application of the formalism displayed in the previous section to the 
particular case known as Taub model. In Bianchi cosmological models |10[ the minisuper- 
space is a three dimensional manifold parametrized by two parameters (/3 + ,/3_) measuring 
the spatial anisotropy and a parameter a measuring the volume of the Universe (Misner 
parametrization). The Hamiltonian constraint for minisuperspace models has the form 

H = e 3a {-p 2 a +pl+ p 2 _ + e~ 4a [V ((3 + , /?_) - 1]} (15) 

where (p a ,p + ,p_) are the momenta canonically conjugate to (a, /3 + ,/3_) and the potential 
V (P + ,{3-) depends upon the particular Bianchi model. 

The Taub model is a particular case of the Bianchi IX for (3_ = 0, p_ = 0. For this case 
the resulting Hamiltonian constraint is 

H=-pl+p%+ 12n 2 e- 4Q (e- 8 ^ - 4e" 2/3 +) (16) 

scaling the Hamiltonian with the factor e~ 3a . If we define the variables u and v by 
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t = Arc sinh j^-^-e- 61 ' ) (19) 
Pt 

36 



a = v - 2u (17) 
/?+ = u - 2v 

the resulting Hamiltonian is (multiplied by 1/6) 

H = I (p 2 v + 36ttV 2 *) -\{pl + 144vr 2 e 6 ") (18) 

By means of this coordinate transformation we could in fact separate in the Hamiltonian 
a subsystem depending on just one pair of canonical variables, which will work as a clock 
for the other subsystem. A global time is defined by means of the canonical transformation 
defined in (0) . In our case the new variables are 

Pv 

6tt 

1 (pi + 36ttV 2 *) (20) 
The generator of this transformation is 

Fi(v,t) = -7re 6u sinh£ (21) 
The Hamiltonian in the new variables results to be 

H = 6 P 2 t ~\{pI + 144vr 2 e 6 ") (22) 
This expression can be factorized in order to obtain a Hamiltonian linear in p t , so giving 
H = (VEpt + ^v^ + ^V 3 ") (V6p t - ^Pl + 1447r 2 e 6 ^ (23) 

The constraint H = is fulfilled if one of the factors is null on the constraint surface. 
The other factor has, on the constraint surface, a definite sign, so playing the role of the 
factor / defined before. The scaled Hamiltonian is 

H 1 



with 



H = -^j=p t + -y/p* + 1447r 2 e6« =p t + h u 



( n: 1 



/= [VQpt-^P 2 u + U^ 2 e &u 
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B. Quantization 



In order to quantize the reduced system we will make the substitution p t 



—i 



of 



Pu - *Q,„ 

equation 



— iS- and impose the constraint H^/ (t, u) = yielding the following Schrodinger 



(25) 



Inserting solutions of the form \I/ (t, u) = (f)(u) e 1 V 6 we obtain a modified Bessel equation 
for the function (f> (u) . The solutions of this equation are the modified Bessel functions 



<P (it) = CK 2iy r £ (4vre 3u ) + DI 2lV ~ £ (4vre 3u 



(26) 



The functions (47re 3lt ) should be discarded because they diverge when u — >• oo (clas- 
sically forbidden zone). The solutions corresponding to the quantization of the reduced 
system are therefore 



(27) 



On the other hand the Wheeler-De Witt equation associated with the Hamiltonian ( |i~8|) is 



6 \ dv' 
whose solutions are 



+ 36?re 



2A2v 



i / d 2 



6 \ du : 



2J6u 



ip (v, u) = 



(28) 



ip (v, u) = [AK iV - £ (ixe^) + BI lV - £ (vre 6 ^)] [CK 2iV ~ £ (4vre 3 ") + DI 2i ^ (4ne 3u 



(29) 



It would not be correct to impose the same kind of boundary conditions used to discard 
the functions I 2 i^ (47re 3u ) in the quantization of the reduced system. The variable v is not 
a dynamical variable but the variable associated with the clock of the system. It is by no 
means obvious that the physical solutions should go to zero in the classical forbidden zone. 

In order to select the physical solutions we will try to apply the "quantum canonical trans- 
formations" defined in (|14"1) to the solutions of the Wheeler-De Witt equation. The physical 
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solutions will be those whose transformed functions are the solutions of the Schrodinger 
equation e™*^'. We will begin by transforming the functions which go to zero in the clas- 
sically forbidden zone, i.e., the functions 0(f) = K i ^( / Ke 6v ). The transformed functions 
are 



+00 



$(t) = N I dvGne*" coshte 

7TN 



iire 6v sinh t 



e 2 e 



-iy/et 



(30) 



4 sinh (^f) cosh(^) 

In this way it is manifest that the transformation of the functions O (v) = K iy ^ (7te 6v ) do 
not give the solutions of the Schrodinger equation e~ % ^ 1 . On the contrary they correspond 
to a combination of positive and negative energies states 



ttN 



4 sinh coshf^# 



e 2 e 



(31) 



By transforming the right side of ( pT| ) one should recover the original function (v) = 
Kiy/e (ne^) , so one obtains the factor N = which does not depend on the energy. 

As the functions (v) = (ire ev ) are not definite energy states, we will apply the 
transformation fll4]) to the other subspace of solutions, i.e., to the functions I±i^ (ne 6v ) , 
which diverge in the classically forbidden zone. The resulting integral has the form 



I f+00 

$(*) = -=/ dv6ire 6v coshte 1 



i7re 6 " sinht 



(32) 



This integral diverges unless one gives an imaginary part r\ to t. Replacing t — > t + i^ in fl32|) 
one can actually perform the integration. Replacing t — ► t — i| in the result of the integral 
one can go back to the original real time variable, obtaining the correspondence 



(33) 



The functions 7 iv ^ (7re 6w ) and (7re 6w ) do represent then the positive and negative energy 
states respectively []. In this way we could establish which subspace of the whole space of 



1 This result can be verified by testing the consistence of ( |33|) with ( pi]) and the expression [11] 
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solutions of the Wheeler-De Witt equation is the physical one. It is remarkable that the 
functions in the selected subspace do not decay in the classically forbidden zone. 



C. Boundary conditions 

It would be interesting if one could define certain boundary conditions which would not 
rely on the fact that one knows how to reduce the system, and which would select the 
same physical functions imposed by the quantization of the reduced system. Let us start by 
considering the following Hamiltonian 

H = p{ + V ( gi ) - h qft PlM ),n = 2,..,N (35) 

and suppose that <p (q 1 , qyj = (qi) <fi (qyj is the solution of the Wheeler-De Witt equation 
associated with this Hamiltonian. Instead of performing a canonical transformation so that 
the new Hamiltonian is quadratic in the new momentum p t , we will study the wave functions 
obtained by solving the Wheeler-De Witt equation associated with the Hamiltonian (^) in 
the region L where V (qi) tends to zero. In this region the Hamiltonian is 

H = p 2 1 -h^(q^p^) (36) 

The solutions of the quantum-mechanical system corresponding to the sheet in which p\ 
is negative (positive energy solutions) will be combinations of tp (qi, q^) = 4> (g M ) e~ % ^ qi . One 
would expect that these solutions do coincide with the asymptotic expressions in the region 
L of the functions y? (qi, g M ) = (qi) (f> (g M ), i.e., it would be necessary that 

Q(qi) -> e"^ 91 (37) 



jf.w = f- ( f'- / ; w (34) 

2 sin [vtt) 

In fact, transforming the right side of (f34[) using p^) one does obtain the right side of (|3l|) , verifying 



in this way the coherence of the found correspondences (31, |33|) between both representations. 
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The definite energy solutions will thus be those functions G (qi) which do behave in the 
asymptotic region like a plane wave, ingoing or outgoing. This criterium relies on the fact 
that in the asymptotic region time is (modulus a sign) the variable q\. 

In order to test this criterium in the case of the Taub model let us see the behavior of 
the Bessel functions in the region where v — > — oo (V (qi) — > 0). The asymptotic expressions 



arc 



h{z) 
KAz) 



Ti 



2 sin (vtt) 



r (-v + 1) 



T(v + 1) 



The asymptotic expression for K v (z) was obtained using the formula 
totic expressions are thus 



i\fe 



K) 



(36tt 2 ) — 



(i2)^r(^ + i) 



(38) 
(39) 

The asymp- 
(40) 



71 



2 sin (vtt) 



(36tt 



^^ 2 



(36tt 2 ) 



^^ 2 



(41) 



(i2p^ r {-iy/i + 1) (\2) lV ~ £ r (i^fe + 1) 

In this way we confirm that the functions (vre 6 *') do correspond to a combination 
of positive and negative energy states. The functions I±i^q (ne 6 "") do correspond to states 
of positive or negative energy respectively. The proposed criterium establishes boundary 
conditions with the definite meaning of selecting the positive energy states. The space of 
solutions of these positive energy states can be endowed with the positive definite Schrodinger 
inner product. It is remarkable that the proposed boundary conditions coincides with the 



criterium proposed by Wald [12,1 



IV. CONCLUSIONS 



In this work it was addressed the question of quantizing minisuperspace models by 
studying the particular case known as Taub Model. The main two problems which arise 
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in the canonical approach are the boundary conditions to be imposed on the solution of 
the Wheeler-De Witt equation and the inner product to be defined in the corresponding 
Hilbert space. In the Taub model it is possible to perform a coordinate transformation in 
order to separate in the Hamiltonian constraint a subsystem depending on just one pair 
of canonical variables. In this new set of variables the Hamiltonian has thus the form 
H = h qi (qi,pi) — h q (q^,,Pfj) where the subsystem h qi will work as the clock of the model. 
Performing a canonical transformation it is possible to transform the subsystem h qi in a free 
system h t = pf . This kind of time variables are known as "extrinsic time" because they are 



associated not only with the coordinates but also with the momenta [p],M|. Extrinsic times 
are specially important in quantum gravity because it is not possible to reduce the system 
by identifying an intrinsic time, i.e., a global time variable in the configuration space, as it 
happens for the relativistic particle. The new Hamiltonian H = pf — h q (q^jPn) can be 
factorized in two disconnected sheets. In order to satisfy the constraint H = it is necessary 
that one of this factors goes to zero on the constraint surface. As the other one has thus a 
definite sign on the constraint surface, it is possible to scale the Hamiltonian constraint in 
order to find a Hamiltonian linear in the new momentum p t . This Hamiltonian can be quan- 
tized by means of an ordinary Schrodinger equation. The canonical transformation used to 
reduce the system satisfies the necessary conditions to define the integral transforms which 
relate the wave functions corresponding to the quantization of both Hamiltonian systems. It 
is thus possible to transform the positive energy solutions of the parabolic Schrodinger equa- 
tion in order to find out those solutions of the hyperbolic Wheeler-De Witt equation which 
are the physical ones. Armed with the knowledge of the physical solutions, we tried define a 
criterium to select those solutions without using the fact that the reduced system is known. 
In order to do that it was studied the asymptotic behavior in the free zone (V (qi) — > 0) of 
the solutions of the Wheeler-De Witt equation. We argue that in that area the time is (minus 
a sign) the variable q%. It is thus necessary that the functions $ (gi) behave as an outgoing 
or ingoing plane wave, i.e., as definite energies states. The wave functions satisfying this 
criterium do coincide with the physical functions selected by reducing the system. 
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